Quasinormal modes and thermodynamic phase transitions 
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It has recently been suggested that scalar, Dirac and Rarita-Schwinger perturbations are related 
to thermodynamic phase transitions of charged (Reissner-Nordstrom) black holes. In this note we 
show that this result is probably a numerical coincidence, and that the conjectured correspondence 
does not straightforwardly generalize to other metrics, such as Kerr or Schwarzschild (anti-)de Sitter. 
Our calculations do not rule out a relation between dynamical and thermodynamical properties of 
black holes, but they suggest that such a relation is non-trivial. 

PACS numbers: 04.70.-s, 04.70.Dy 



I. INTRODUCTION 

Kerr-Newman black holes (BHs) in (anti-)de Sitter 
space are characterized by their mass M, dimensionless 
rotation parameter a/M, dimensionless charge Q/M and 
possibly a non-zero cosmological constant A (here and in 
the following we adopt geometrical units, G = c = 1, and 
we write Einstein's equations as G M „ — 3A<? M „.) 

Davies 0, Q pointed out that BHs can undergo a 
second-order phase transition in which their specific heat 
changes sign as a/M and/or Q/M are increased. For 
our considerations below, we only need to recall that the 
phase transition occurs when 
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(1) 



for Reissner-Nordstrom (anti)-de Sitter BHs (a = 0) and 



\MJ 



a/5-1 



(2) 



for Kerr BHs (A = Q = 0). 

BH perturbations due to external fields of different 
spin have been studied for decades 0. After a transient 
phase, the decay of these perturbations can be described 
as a superposition of damped exponentials of the form 
exp(iwi). For a given spacetime (i.e. for a given cos- 
mological constant A) the quasinormal mode (QNM) fre- 
quencies u) = LUR+iujj depend only on the BH parameters 
M , a and Q. They are complex numbers because oscilla- 
tions are damped away by the emission of gravitational 
radiation, and they are usually labeled by three integers: 
uj = uji mn (see [H for a review, and [j| for a summary 
of the properties of the QNM spectrum.) Two integers 
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(I, m) correspond to the "angular quantum numbers" 
of the spin-weighted spheroidal harmonics used to sep- 
arate the angular dependence of the perturbations. For 
spherically symmetric backgrounds, such as the Reissner- 
Nordstrom metric (a = A = 0), perturbations are de- 
generate with respect to the azimuthal index m, but 
this degeneracy is broken when a ^ 0. A third integer 
n(= 0,1,...), called the "overtone number", sorts the 
frequencies by the magnitude of their imaginary parts. 

For Reissner-Nordstrom BHs and n < n c [n c being 
some critical value that depends on I and on the spin s 
of the perturbing field), ujr and wj are usually monotonic 
functions of the charge Q/M. For n = n c the real part 
of the QNM frequency has an extremum as a function of 
the charge. When n > n c , both cor and wrusually be- 
come oscillatory functions of Q/M jy, @, H, Q. A similar 
behavior has been observed for m = QNM frequencies 
of Kerr BHs as function of a/M 0, M, E| (see also [H 
for a generalization to Kerr-Newman BHs.) 

Jing and Pan [l3| recently computed QNM frequencies 
for scalar (s = 0), Dirac (s — 1/2) and Rarita-Schwinger 
(s = 3/2) perturbations of Reissner-Nordstrom BHs. 
They noticed that the first maximum of Mujr(Q/M), 
as computed numerically, matches within ~ 2.5% the 
value of Q/M = ^/i/l ~ 0.866 predicted by Eq. Q 
for Davies' second-order phase transition. They conjec- 
tured that this agreement implies a connection between 
dynamical and thermodynamic properties of Reissner- 
Nordstrom BHs. In this note we show that the observed 
agreement is probably a numerical coincidence, and that 
such a connection (if it exists) should be non-trivial. 

Our results can be summarized as follows: 

(1) The approximate agreement found in Ref. [l3| is usu- 
ally observed for n = n c + 1, not for n = n c . We see no 
compelling reason why this mode should be singled out 
from the QNM spectrum as especially relevant. 

(2) For n = n c + 1, the numerical agreement with Davies' 
thermodynamic phase transition point is significantly 
worse for "electromagnetic-type" or "gravitational-type" 
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perturbations of Reissner-Nordstrom BHs, which were 
not considered in Ref. (l3| . 

(3) Finally, and perhaps more convincingly, we show that 
the conjecture does not hold for (i) integer-spin, m = 
perturbations of Kerr BHs, and (ii) integer-spin pertur- 
bations of Schwarzschild anti-de Sitter (SAdS) BHs. 

In Section [IT] we present details of our calculations. In 
Section IIIII we argue that, if a relation between classi- 
cal BH oscillations and their thermodynamic properties 
exists, it is not as simple as proposed in Ref. [13]. 

II. RESULTS 

A. Reissner-Nordstrom black holes 

Ref. [HI considered scalar (s = 0), Dirac (s = 1/2) and 
Rarita-Schwinger (s = 3/2) perturbations of a Reissner- 
Nordstrom BH. The results presented in their Figure 1 
for s = agree with those previously published in [l2| . 
From the calculations shown in Figure 1 of [13] (or Figure 
1 of [lH, where modes are counted starting from n = 1 
rather than from n = 0) one deduces that, for scalar 
perturbations with I — 0, Mlor has the first extremum 
when n = n c = and Q/M = 0.962. Similarly, for scalar 
perturbations with 1 = 1, the first extremum occurs when 
n = n c = 1 and Q/M = 0.960. Both values are quite far 
from Davies' phase transition point Q/M ~ 0.866. 

TABLE I: "Critical" charge Q/M corresponding to the first 
extremum of Mur(Q / M) for Reissner-Nordstrom BHs. We 
consider only the dominant multipolar component of the ra- 
diation, i.e. I = \s\. For each spin s, in the first row we list the 
critical charge for n = n c (the overtone number for which an 
extremum in Mlor first appears); in the second row, we show 
the charge corresponding to the first extremum for n = n c + l. 



Q/M 


Perturbation 


/ = 


I = 1 


I = 2 


Scalar 


0.962 (0) 
0.879 (1) 






Elect romagnet ic-ty p e 




0.962(2) 
0.940(3) 




Gravitational- type 






0.958(1) 
0.910(2) 



For s = 0, the maximum in Mlur(Q /M) gets in better 
agreement with the phase transition point if we consider 
the next overtone n = n c + 1, as we show by an explicit 
calculation in Table [TJ For the s — I = 0, n = 1 mode, 
and accounting for their different choice of units, Ref. (l3| 
determines the "critical" charge to be Q/M — 0.876. 
This is in good (if not perfect) agreement with our value 
Q/M = 0.879, which has a 1.5% disagreement with the 
phase transition point predicted by Davies. 

Ref. [HI did not consider gravitational (s = 2) 
and electromagnetic (s = 1) perturbations of Reissner- 
Nordstrom BHs. These perturbations are coupled to each 
other by the BH charge, but the scattering problem can 



still be reduced to a pair of wave equations (see discus- 
sions in [3, E3, EH-) The two equations are said to de- 
scribe "electromagnetic-type" and "gravitational-type" 
perturbations, depending on whether they reduce to pure 
electromagnetic or pure gravitational perturbations of a 
Schwarzschild BH in the limit Q/M -> 0. In Table [J 
we list the critical charge for the dominant multipoles of 
"electromagnetic type" and "gravitational type" pertur- 
bations. From the tabulated values it is quite clear that 
the agreement with Davies' phase transition point does 
not get better even if, following Ref. [H| , we consider the 
QNM with n = n c + 1 rather than that with n = n c . 

In our opinion, this is evidence that the conjectured 
correspondence is only a numerical coincidence. The 
value of the charge for which Mlor first has an extremum 
is around Q/M ~ 0.96 for all integer spins s = 0,1, 2, and 
this is quite far from the thermodynamic phase transi- 
tion point Q/M = 0.866. One could still argue that the 
dynamical and thermodynamic properties of BHs are not 
simplest when one considers coupled electromagnetic and 
gravitational perturbations: see e.g. [3] for a discussion. 
To show that this is not the only reason for the discrep- 
ancy with Davies' critical point, in the remainder of this 
paper we show that the conjectured correspondence does 
not hold for other metrics as well. In particular, below 
we consider Kerr and SAdS BHs. 



B. Kerr black holes 

For Kerr BHs, according to Eq. |[5J) the thermody- 
namic phase transition corresponds to a rotation param- 
eter a/M ~ 0.786. For moderate values of n, only Kerr 
perturbations with m = are oscillatory functions of 
a/M 0, fTTj] . For this reason, in our discussion we only 
consider QNM frequencies with m = 0. 

TABLE II: Dimensionless angular momentum parameter a/M 
corresponding to the first extremum of Mua(a/ M) for Kerr 
perturbations with m = 0. We only consider the two lowest 
allowed values of I: J = |s|,|s| + 1. In parentheses we list n c , 
the overtone number for which an extremum first appears. 



a/M 


Perturbation 


1=0 1=1 1=2 


I = 3 


Scalar 


0.819 (0) 0.898(1) 




Electromagnetic 


0.922(1) 0.925(2) 




Gravitational 


0.947(2) 


0.937(2) 



The first occurrence of local extrema in Mu>n(a/M) 
for m = and different values of I and s is listed in Table 
HIl For n — n c there is no agreement with the ther- 
modynamic phase transition point predicted by Davies. 
As in the Reissner-Nordstrom case, we verified that the 
agreement would not improve if we considered the next 
overtone n = n c + 1. For example, for s = I = m = and 
n = 1 the first extremum in Mojn(a/M) occurs when 
a/M = 0.503. 
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C. Schwarzschild-(anti-)de Sitter black holes 

Setting Q = in Eq. fl) wc get = (9/4) 2 AM 2 + 3/4, 
r + = 3M/2. A phase transition occurs when a/|A| M = 

-2/(3 V3) and v / lA|'"+ = VV3 ~ 0.577 for negative 
cosmological constant (anti-de Sitter space.) For positive 
A (de Sitter space), there is no phase transition. 



TABLE III: BH horizon radius corresponding to the first 
extremum of u>n(r+)/ i/\A\ for SAdS BHs. The number in 
parentheses is n c . 



s/\Mr+ 


Perturbation 


1 = 1 = 1 1 = 2 1 = 3 


Scalar 


0.39 (0) 0.61(0) 0.66 (0) 0.78 (0) 


Electromagnetic 


0.92 (0) 0.92 (0) 


Odd gravitational 


0.75 (1) 1.26 (1) 


Even gravitational 


0.73 (0) 1.04 (0) 



The perturbation equations for SAdS spacetimes have 
been worked out in the literature, and QNM fre quen cies 
have been studied in many papers (see e.g. (l5l. ll6L Il7l 
fbsl].) We used numerical codes we developed in the past 
to look for extrema in the QNM frequencies. We have 
two choices to make the frequency dimensionless: 

(i) We can rescale by the BH mass and compute 
Mu>n(r+). This function has no extrema for SAdS BHs. 

(ii) We can rescale frequencies by the cosmological con- 
stant A. Then u)n(r + )/ \/fA| does have an extremum for 
the values or r+ listed in Table HTT1 compare with Fig. 1 
of Ref. [IT] . It is clear from the Table that these extrema 
are strongly dependent on I and s, and that in general 
they do not agree with Davies' phase transition point. 
The fact that QNM frequencies do have an extremum 
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III. CONCLUSIONS 

In the last few years, remarkable relations between 
classical and thermodynamic properties of black objects 
have been uncovered. For instance, a correspondence 
between classical and thermodynamic instabilities of a 
large number of black branes conjectured by Gubser and 
Mitra [H, HJ was proved by Reall [H| (see [H for a 
review.) Manifestations of this duality are expected to 
appear in the QNM spectra. Indeed, some indications 
that phase transitions do show up in the QNM spec- 
trum were provided in specific cases by various authors 
[H 0, [H, HE H3- The main purpose of this note is 
to show that the correspondence is not as simple as pro- 
posed in [l~3l |. 
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